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Pauli paramagnetic susceptibility x is unaffected by the electron-phonon interaction in the Migdal- 
Eliashberg context. Fullerene compounds however do not fulfill the adiabatic assumption of Migdal's 
theorem and nonadiabatic effects are expected to be relevant in these materials. In this paper we 
investigate the Pauli spin susceptibility in nonadiabatic regime by following a conserving approach 
based on Ward's identity. We find that a sizable renormalization of \ due to electron-phonon cou- 
pling appears when nonadiabatic effects are taken into account. The intrinsic dependence of \ on the 
electron-phonon interaction gives rise to a finite and negative isotope effect which could be experi- 
mentally detected in fullerides. In addition, we find an enhancement of the spin susceptibility with 
temperature increasing, in agreement with the temperature dependence of \ observed in fullerene 
compounds. The role of electronic correlation is also discussed. 

PACS number(s): 74.70.Wz, 7I.38.-k, 63.20.Kr 



I. INTRODUCTION 

The Pauli susceptibility, x? of wide-band metals is usu- 
ally considered as unrenormalized by electron-phonon eL 
fects. As pointed out some time ago by Fay and AppelJil 
X is actually affected by the electron-phonon coupling, 
the effect having magnitude of order Xujp^/Ep, where A 
is the electron-phonon coupling, u p h is a typical phonon 
frequency and Ep is the Fermi energy. For wide-band 
metals, w p h <C Ep and the electron-phonon renormaliza- 
tion of x is therefore negligible. This result basically fol- 
lows from the observation that the lowest order electron- 
phonon correction to \ i^a vertex diagram which, accord- 
ing to Migdal's theoremjj is of order \uj p h/ Ep. Hence the 
absence of phonon effects in \ is j us t a manifestation of 
the electron-phonon adiabaticity of normal metals. 

This situation may be somewhat, different for weak 
ferromagnetic metals like ZrZn2.EHj In this case, the 
electron-phonon corrections compete with the exchange 
term leading to a possible phonon-dependent ferromag- 
netic transition temperature. Experiments have how- 
ever reported only small isotope effects in ZrZn2 with 
quite large error barsjj leaving the problem of phonon- 
corrected weak-ferromagnetism essentially unresolved. 
The recent discovery of giant isotope shifts in manganites 
proves however the existence of ferromagnetic materials 
with important electron-phonon effects in their magnetic 
properties a 

The possibility of having electron-phonon dependent 
spin susceptibility is however not correlated exclusively 
to the vicinity of magnetic instabilities. In fact as long 
as Xujph/Ep is not negligible, as in very narrow band 
systems, the electron-phonon vertex contributions are no 
longer unimportant so that x is expected to acquire a 
phonon renormalization. The most promising candidates 
for the observation of this effect are the Ceo-based ma- 



terials. Fullerene compounds have in fact phonon modes 
extending up to 0.2 eV and Fermi energies of about 0.3 
eV.El Therefore, ujph/Ep is in principle large suggesting 
that Cgo-based metals could be labeled as nonadiabatic 
systems. 

Recently, the finding of a superconducting transition 
at T c — 52 K in hole doped Ceo single-jCrystals has rised 
a renewed interest in these materials!] Nonzero isotope 
effects and other properties strongly indicates that su- 
perconductivity in fullerides is driven by electron-phonon 
interaction.^ However, the description of superconductiv- 
ity of RbaCgo within the traditional Migdal-Eliashberg 
(ME) theory is found to be inconsistent with respect to 
the adiabatic hypothesis Xuip\JEp <C 1 which is at the 
basis of the ME theory itself.H Instead, by relaxing the 
adiabatic hypothesis, a generalized formulation which in- 
cludes nonadiabatic channels in the electron-phonon in- 
teraction provides a more self-consistent decription of su- 
perconductivity, and suggests that the key ingredient for 
the high values of T c in Cgo-based materials is a con- 
structive nonadiabatic interference rather than strong 
electron-phonon couplings n a 

In principle, the hypothesis that superconductivity in 
fullerene compounds is enhanced essentially by nonadi- 
abatic electron-phonon effects can be sustained by the 
observation of independent signatures of nonadiabaticity. 
In this respect and according to what we have pointed 
out at the beginning, Pauli susceptibility is a quantity 
where such signatures could be found. In this paper, we 
provide extensive calculations of \ beyond the adiabatic 
limit by including nonadiabatic effects at different stages 
of a perturbation theory in Xojp^/Ep. We find that when 
w p h/ Ep is no longer negligible, (i) the paramagnetic spin 
susceptibility can be considerably reduced with respect 
to the adiabatic limit, (ii) it acquires a negative isotope 
effect and (iii) a possible anomalous temperature depen- 
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dence at constant sample volume. These features are sig- 
natures of nonadiabatic electron-phonon interaction and 
prediction (ii) and (iii) are susceptible of experimental 
verification. 



II. PAULI SUSCEPTIBILITY BY WARD'S 
IDENTITY 

As remarked in the introduction, the effect of electron- 
phonon interaction on spin susceptibility has already at- 
tracted some interest in the past mainly in relation to 
weak antiferromagnets-p. However, different approaches 
led to different results,Em reflecting the lacking of a con- 
trolled theory. 

In our paper we use the functional formalism based 
on the Baym-Kadanoff technique to derive a conserving 
derivation of Pauli susceptibility valid for both electron- 
phonon and electron-electron interactions. The Pauli 
susceptibility \ is calculated by the knowledge of the spin 
vertex function which is related to the self-energy via a 
Ward's identity. 

Following the Baym-Kadanoff technique, we introduce 
an external magnetic field H coupled with the elec- 
trons which induces a magnetization M . The interac- 
tion hamiltonian describing the coupling of H with the 
electron spins is: 



Hi 



k,<r 



(1) 



where h = n^H and cj. (ck CT ) are creation (annihilation) 
fermionic operators for electrons with wave number k and 
spin index a = ±1. The electron magnetization M due to 

(|l|) is given by M = [i B J2a an °- where n a = Sk( c L c kcr) 
and (■ • •} denotes the statistical average. We express now 
the magnetization M in terms of the finite temperature 
single electron propagator 



(2) 



where T T is the time ordering operator and r is the imag- 
inary time. Since (cj^Cko-) = G a (k, 0~), the magnetiza- 
tion can be expressed as: 



M = ^B^^(k,0-) 



k.tT 



n k,(T 



(3) 



In the above expression, G a (k, n) is the thermal electron 
propagator which satisfies the following Dyson equation: 

G-^n) =ioj n -e(^)+fi + ha-E a (k,n), (4) 

where uj n = (2n+l )ttT, n = 0, ±1, ±2, . . ., are Matsubara 
frequencies, e(k) the electron dispersion, /i the chemical 
potential and S (T (k, n) is the electronic self-energy due 



to the coupling to phonons and to the electron-electron 
interaction. The spin susceptibility \ is formally given 
by M — xH where \ = [dM/dH]o is the derivative of 
the magnetization at zero field. Hence, from Eq. (|J), a 
general expression for \ is the following: 



dG a (k, n) 
dh 



= -2^T^]TG(k,n) 2 r(k,n) 



(5) 



where G(k, n) is the electron propagator for zero mag- 
netic field which satisfies the H — > limit of Eq. fl) : 



G _1 (k, n) 



e(k)+/i-E(k,n). 



(0) 



In the second term of Eq. (||) we have introduced the 
spin- vertex function: 



r(k,n) = -5> 



dG-^Kn) 
dh 



(7) 



In equations (||) and (^), the notation [• • -]o indicates 
that the quantity in brackets must be calculated for zero 
magnetic field. Plugging Eq. (|J) into (0) the resulting 
vertex T satisfies the Ward's identity: 



r(k,n) = i-~5> 



dS CT (k, n) 



dh 



(8) 



and, since E CT (k, n) is a functional of the electron prop- 
agator, a self-consistent relation between (^) and (||) is 
obtained which permits to calculate the spin suscepti- 
bility. At this point, the spin susceptibility \ can be 
calculated once the electron self-energy and its magnetic 
field dependence is known. 



III. NONADIABATIC PAULI SUSCEPTIBILITY 

Electron-phonon interaction is usually neglected in the 
calculations of the Pauli susceptibility although it can 
strongly renormalize other physical quantities .113 Indeed 
it can be shown that the electron-phonon self-energy de- 
pends on the external magnetic field h at least as 



lim E cp (k,n) 



hO(uj ph /E F ), 



(9) 



where w p h characterizes the phonon frequency scale and 
E-p is the Fermi energy. In wide-band materials the adia- 
batic ratio m-p^/E-p <C 1 and the electron-phonon contri- 
bution to x can be consequently disregarded. From Eq. 
(H) , we see that electron-phonon effects appear in \ onr y 
at a nonadiabatic level. Hence, any evidence of electron- 
phonon effects on the Pauli susceptibility is thus a direct 
proof of a nonadiabatic electron-phonon coupling. 
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FIG. 1. Electronic self-energy within the non-crossing ap- 
proximation. 

In order to properly include electron-phonon interac- 
tions in x we need then to explicitely specify the func- 
tional form of the nonadiabatic electron-phonon self- 
energy. Nonadiabatic effects enter in a twofold way in the 
electron-phonon self-energy: finite bandwidth effects and 
vertex diagrams. These two kinds of effects, of course, are 
of the same order and there is no justification for neglect- 
ing vertex corrections with respect to finite bandwidth 
effects. 

In this paper we consider two approximation schemes 
for the electron-phonon self-energy. The first one is 
essentially the mean-field theory which corresponds to 
the non-crossing approximation. It is diagrammatically 
equivalent to the Migdal-Eliashberg electron-phonon self- 
energy without however assuming an infinite electron 
bandwidth compared to the relevant phonon energies. At 
this level only finite bandwidth nonadiabatic effects are 
considered. The second one includes first order electron- 
phonon vertex corrections as well as finite bandwidth ef- 
fects inpthe framework on the nonadiabatic Fermi liquid 
pictured Both approximation schemes reduce to the adi- 
abatic Migdal-Eliashberg limit for w ph /_E F — » 0. 



A. Non-crossing self-energy 

The self-consistent non-crossing approximation ne- 
glects the vertex corrections in the electron-phonon self- 
energy and for utp^/Ep <C 1 reduces to the ME theory of 
the electron-phonon coupled system. The diagrammatic 
representation of the electron self-energy is shown in Fig. 
[l] where the wiggled line represents the phonon propaga- 
tor and the dashed line the electron-electron Coulomb 
repulsion. The corresponding compact expression of the 
non-crossing self-energy reads: 



M*) = J2 [ V( - k - fc ') + Ie- 1 ^] G a {k' 



(10) 



where k and k' are fermionic four-vectors defined as k = 
(k, iw n ) and k' = (k',iu> m ). Moreover J2k = ~ T J2nJ2k 
and V(k-k') = \g(k-k')\ 2 D(k-k'), where ff(k-k') is 
the electron-phonon matrix element and D(k — k') is the 
phonon propagator. / is exchange Coulomb interaction 
which gives rise to the Stoncr enhancement factor. 

By introducing Eq. ([l0]) into the expression of the 
spin- vertex function (g) and using Eq. (f?) we obtain: 

dG a (k') 



T(k) = 1-J2 [V(k - k') + Ie-**»°-] J2 \ 

k' a 

= 1 + - k ') + /]G(fc') 2 r(fc')- 



dh 



o 

(11) 
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FIG. 2. Diagramatic self-consistent expression of the spin 
vertex in non-crossing approximation. 

The non-crossing approximation for the self-energy leads 
therefore to a self-consistent ladder equation for T(k, n)EJ 
(Fig. ||). If we consider dispersionless phonons of fre- 
quency uiq, the above ladder equations can be rewritten 
in the extended notation as: 

T(k, n) = 1 - T J2 [ls(k - k r )\ 2 D(n - m) + i] 

k' ,m 

xG(k',m) 2 r(k',m), (12) 

where 



D(n — m) — 



(iuj m - iuj m ) 2 - u 2 



(13) 



The electron propagator appearing in Eq. (|llj) satisfies 
the Dyson equation (Q) with the zero field limit of the 
self-energy (|i~0|): 

E(k,n) = -T^2 |.g(k-k')| 2 J D(n-m)G(k',m) 

k' ,m 

+IT G ( k '> m)e-^°~ , (14) 

k' ,m 

In the simplest case of k-independent Coulomb repul- 
sion here considered the second term on the right side 
of Eq. (|l4|) gives rise just to a constant term which can 
be absorbed into a redefinition of the chemical potential 
pi — > p. We can then neglect it since we shall consider 
only half-filled systems for which we set pi = 0. 

In an isotropic system the angular dependence of the 
self-energy and of the vertex function in equations ( |T^ ) 
and (|lj) is negligible and it can be dropped. Therefore, 
following the same procedure reported in Ref. O, we re- 
place the electron-phonon interaction |g(k — k)P by its 
average over the Fermi surface: 



|. 9 (k-k')| 2 -<<|. 9 (k-k')| 2 )>FS^ 2 , (15) 
£ k . k , | 5 (k-k')| 2 <5[e(k)M e (k')] 



where 
((lff(k-k')| 2 )) FS 



Ek,k'^(k)Me(k')] 



(16) 



In this way, the electron self-energy becomes indepen- 
dent of the momentum, S(k, n) = £(n), and the electron 
propagator becomes at half-filling: 



G(k,n) = 



1 



W n - e(k) • 



(17) 
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FIG. 3. Vertex corrected electron-phonon self-energy. 



where we have set S(n) = iuj n — iW n - By using a con- 
stant density of states, Nq, the momentum summation 
in Eq. ( |l4| ) is transformed as — > de and, 

under integration over e, the renormalized frequency W n 
becomes: 



W n = uj r , 



2 / .Bp \ 
XttT D(n — m) — arctan f I 



(18) 



where A = Nog 2 is the electron-phonon coupling con- 
stant. By using ( [l5|) and (|l7|), also the spin- vertex func- 
tion ( |l2| ) becomes momentum independent and, by fol- 
lowing the same steps as above, it reduces to: 



r(n) = 1 + [\D(n -m) + I] 



2E F 



Wl+E 2 F 



r(m). (19) 



Finally, since the self-energy and the spin- vertex function 
depend only on the frequency, the summation over k can 
be readily performed in (||) leading to: 



X(T) = XpTY, w 2 2 ^ F F 2 r W> 



W 2 + E 2 



(20) 



where \p = 2^|Ao. The spin susceptibility is then ob- 
tained by the solution of equations ( |l8| ) , ( [l9| ) , and (^) . 



B. Vertex corrected self-energy 

In the vertex correction approximation, the electron- 
phonon self-energy E cp is modified with respect to the 
ME one by the inclusion of the first electron-phonon ver- 
tex diagram as shown in Fig. ^ By making use of the 
condensed notation introduced in Sec. [II A, the vertex 



corrected self-energy can therefore be written as follows: 



l+J2v(k-q)GAq-k + k')G a (q) 



(21) 



where q = (q, iu>i). The derivative of £° P (A;) with respect 
to h = fisH calculated at zero magnetic field is: 



dK p (k) 
dh 



dG a (k') 



dh 



v ( k - k')V(k - q)G(q -k + k')G(q) 

k' ,q 



dG a {k') 
dh 
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FIG. 4. Diagramatic self-consistent expression of the spin 
vertex in the vertex corrected theory. 



J2V(k-k')V(k-q)G(k')G(q) 



k' : q 



dG a (q-k + k') 



dh 



J2 V ( k ~ k')V(k - q)G{k')G{q -k + k') 

k' ,q 



dG a (q) 
dh 



o 

(22) 



and by re-arranging the 4-vectors indeces and using (Q) 
and (Mj the spin-vertex equation reduces to: 

T(k) = l + J2 [^(*> fc ') + J ] G(k') 2 T(k r ), (23) 



where the electron-phonon kernel V(k, k') is: 



V{k,k') = V(k-k') 



1 + V(k - q)G(q)G{q -k + k') 



V(k - q)V{q - k')G(q)G{k + k' - q). (24) 



in terms of 



A graphical representation of ( |23| ) and (£ 
Fcynman diagrams is shown in Fig. 4. 

In comparison with Fig. ^, the kernel of the spin sus- 
ceptibility resulting from the vertex corrected self-energy 
is modified by additional electron-phonon contributions 
represented by vertex and cross diagrams. This set of di- 
agrams is therefore quite similar to those encountered in 
the theory of nonadiabatic superconductivity. They how- 
ever differ in the orientation of the electron propagators 
(particle-hole rather than particle-particle contributions) 
leading to a different momentum dependence of the cross 
term. The evaluation of the kernel (|2^) can therefore 
be carried on by following a scheme. .similar to the one 
already employed in previous works Jj'Ej 

Let us first evaluate the zero-field self-energy entering 
the electron propagators in Eq. (p4|). As discussed be- 
fore, the Coulomb self-energy provides just a constant 
term which can be absorbed into the definition of the 
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chemical potential. The total self-energy reduces there- 
fore to the electron-phonon one which from Eq. (|2l]) can 
be written in the compact form: 

E(k,n)= -T^^|.g(k-k')| 2 J D(n-m) 

m k' 

x [1 + P(k, k'; n, to)] G(k', m), (25) 
where we have introduced the vertex function given by: 

P(k, k'; n, to) = -T^ £ | 5 (k - q)| 2 #(n - Z) 
I q 

xG(q-k + k',Z-n + m)G(q,Z). (26) 

In similar way as in non-crossing approximation, we 
eliminate the angular dependence of the self-energy by re- 
placing the whole electron-phonon matrix element, which 
takes into account the vertex correction, by its average 
over the Fermi surface: 

|.g(k-k')| 2 [l + P(k,k';n,m)] _ 

((| 5 (k-k')| 2 [l + P(k,k';n,m)])) FS 

= g 2 [l + XP(Q c ;n,m)} (27) 

The analytic expression of P(Q c ;n,m) and its explicit 
derivation is provided in Appendix 0. The parameter Q c 
is a dimensionless momentum describing the relevant mo- 
mentum scale of the electron-phonon scattering process. 
In conventional metals, as the low temperature supercon- 
ductors, the maximum exchanged phonon momentum q c 
is about the Debye vector q^> and Q c ~ qx>/2kF ~ 1- 
It is well known, however, that in systems with low 
density of charge carriers, as fullerene's compounds and 
cuprates, the electrons scatter only with charge fluctu- 
ations of long wavelength because of the strong elec- 
tronic correlation. In fact, in strongly correlated systems, 
the electrons are surrounded by giant correlation holes 
which sujpHress charge density fluctuations with large 
moment a. The relevant exchanged phonon scale is 
thus quite smaller than the Debye vector: q c < q^ 
and Q c — q c /2kp < 1. As a consequence the effective 
electron-phonon interaction is suppressed for momenta 
larger than some cut-off q c which depends on the filling 
factor and on the Coulomb repulsion. We modelize this 
situation by introducing a sharp momentum cut-off q c . 
It can be used as a free parameter to distinguish between 
weak and strong correlation cases, where q c is small for 
strong correlated systems and of order one for weakly 
interacting electrons. 

As usual, the self-energy effects can be expressed in a 
convenient form by means of the renormalized frequencies 
W n = fc> w [l — S(n)/(ztt>„)] obtained by integrating Eq. 
(|25|) over the energy. In the vertex corrected theory they 
fulfil the self-consistent relation: 

W n = w„ - ttT c ^ A[l + XP(Q C ; n, m)]D{n - to) 



By using the momentum independent self-energy fl28|), 
the spin-vertex function (|23|) can be rewritten as: 



r(k',m) 



r(k, n) = i-rJ2 I 1 + V Q*, k '; ^ m ) 



k'.m 



6(k')] 2 ' 

(29) 



where the nonadiabatic electron-phonon spin kernel 
V(k, k'; n, to) is given by: 

F(k,k';n,m) = \g(k - k')\ 2 D(n - m) 

x [1 + 2P(k, k'; n, to)] + G(k, k'; n, to). 

(30) 

Here P(k, k';n, m) is again the first vertex correction 
given in Eq. (p6|) and G(k, k';n, m) is the cross cor- 
rection whose the explicit expression is given below: 

G(k, k'; n,m)=T^2Yl l^( k - <l)\ 2 \9(<l - k ')| 2 
i q 
D(n - l)D(l - m) 



\iWi - e(q)]\iW n 



+m— Z 



e(k + k' - q)] ' 



(31) 



Coherently with the approximations performed on the 
self-energy and with the Ward's relation Eq. (||), we 
evaluate T(k, to) by replacing the kernel ( |30| ) by its mo- 
mentum average over the Fermi surface: 

F(k,k'; n,m) -> {{V(k, k'; n, m)» PS = V{Q c ;n,m), 

(32) 

where 

V (Q c ; n, to) = XD(n - m) [1 + 2\P{Q C ; n, to)] 

+ X 2 C(Q c ;n,m). (33) 

The explicit expression of C(Q C ; n, to) can be also found 
in Appendix |A[ 

The final expression of the ladder vertex equation be- 
yond the adiabatic approximation is readily obtained 
from Eqs. (p9j)-(|33|). The result of the integration over 
the electron energy gives: 



r(n) = 1 - XT^2 l 1 + v (Qc, n, to) 



2E f 



Wl + E 2 F 



T(m). 



(34) 



Finally, the spin-susceptibility in the vertex corrected ap- 
proximation is obtained by plugging Eq. (|34|) into Eq. 
©■ 



IV. RESULTS 



2 (E F 
x — arctan 

7T \W m 



(28) 



We are now in the position to calculate the Pauli sus- 
ceptibility x an d to evaluate the effects on x °f t ne 
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FIG. 5. Spin susceptibility as function of the adiabatic pa- 
rameter uo/Ef and of the electron-phonon coupling A for 
No I = 0.4. Dashed lines represent the spin susceptibility 
for non-crossing approximation, solid lines for the vertex cor- 
rected theory (from the lower to the upper line: Q c = 0.1, 
Q c = 0.3, Q c = 0.5, Q c = 0.7, Q c = 1.0). 



electron-phonon interaction, both in the non-crossing 
approximation and in the vertex corrected theory. Of 
course, when the adiabatic parameter ujq/Ep or the 
electron-phonon coupling constant A are turned to zero 
X would reduce to the simple Stoner enhanced suscepti- 
bility: 

lim v(T)= lim v(T) = X ° ( ^ ( ) , (35) 

where Xo(T) is the free-electron Pauli susceptibility: 

1 



CT) = -2/4T£Eli 

n k 1 

= Xp[1- 2/(^)1, 



s(k)]' 



(36) 



where J(Ef) is the Fermi distribution function at Ef- 

In the following we will denote the zero electron- 
phonon limit in Eq. ( |35| ) as x° c and it would be used 
as comparison element to evaluate the effects of the 
electron-phonon interaction. 

In Fig. [|we plot the total spin susceptibility (electron- 
electron + electron-phonon scattering) as function of the 
electron-phonon coupling and of the adiabatic parameter 
for zero temperature. Dashed lines are the results ob- 
tained within the non-crossing approximation while the 
solid lines are the data for the vertex corrected theory. 
For this latter case we show the results for different val- 
ues of the momentum cut-off Q c (Q c — 0.1, 0.3, 0.5, 0.7, 
1.0). 

The first main result of Fig. || is that the inclusion of 
the electron-phonon coupling, in the nonadiabatic regime 
loq/Ef > yields a sensible reduction of x with respect 
to the pure electronic spin susceptibility. As expected 
this effect vanishes as A — > (right panel) or ojo / Ep — » 
(left panel). Note that both the non-crossing and vertex 
corrected theories yield similar reduction. This is quite 
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FIG. 6. Isotope effect on the spin susceptibility as func- 
tion of loo/Ef (A = 0.7, NqI = 0.4) and as function of A 
(uo/E-p = 0.7, Nol = 0.4). Solid lines and dashed line as in 
previous captions. 



different from the situatkaipCricountered in the supercon- 
ducting pairing channel Jajj'E3 where the effect of the ver- 
tex corrections is much stronger and highly dependent on 

Qc- 

The results of Fig. || suggest that some care is needed 
in estimating the value of the bare density of states from 
paramagnetic susceptibility measurements. Indeed, our 
analysis shows that x is n °t simply related to the density 
of states by Eq. (|3^) . Namely, disregarding the electron- 
phonon effects would lead to a substantial understima- 
tion of the bare density of states from a spin susceptibility 
measurement as long as the electron-phonon interaction 
is in the nonadiabatic regime. 

A more clear signature of the nonadiabatic effects is 
provided by the isotope dependence of the spin suscep- 
tibility. In Fig. H we report the numerical calculations 
of the isotope coefficient a x — — dlogx/dlogM, where 
M is the ion mass, as a function of the adiabatic ratio 
uJo/Ep and of A. Both the non-crossing (dashed lines) 
and the vertex corrected (solid lines) theories predict neg- 
ative values of a x . Compared to the non-crossing data, 
the vertex corrected results show for small values of Q c 
a stronger dependence on lu /Ef and A leading to a x of 
about —0.1. The observation of the isotope effect, which 
should be absent in metals fulfilling the ME framework, 
could be therefore a stringent evidence of a nonadiabatic 
electron-phonon coupling. Note that a previous analysis 
of experimental data of some superconducting proper- 
ties has demonstrated the failure of the ME theory for 
RbaCeo, pointing out the breakdown of Migdal's theo- 
rem in fullerides. In this respect, a measurement of a x in 
Ceo compounds would represent a direct and independent 
test of the relevance of nonadiabatic electron-phonon in- 
teraction. 

From a qualitative point of view, the zero tempera- 
ture behaviors of both the non-crossing and vertex cor- 
rected theories give substantially similar results. A in- 
teresting differentiation among the two approximations 
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FIG. 7. Temperature dependence of the spin susceptibility 
for different values of ujo/E-p and A and Nol = 0.4. Solid lines 
and dashed line as in previous captions. 



arises in the temperature dependence of x- I n Fig- @ 
we report the calculated temperature dependence of the 
Pauli susceptibility for different values of aj Q /Ep and A. 
The non-crossing approximation (dashed lines) yields a 
monotone decreasing of \ as the temperature T increases. 
This is basically due to the phonon cut-off in the lad- 
der equation for the susceptibility, Eq. (|l9|). On the 
other hand, the vertex corrected \ (solid lines) has a 
richer temperature dependence which becomes more im- 
portant as XuJo/Ep increases. Starting from T = 0, the 
basic feature is represented by an initial increase of x 
with T followed by a decreasing for larger temperatures. 
Although the decreasing part is rather Q c independent, 
the initial increase of x is steeper for lower values of the 
momentum cutoff. Note for example that for A = 0.7, 
loo/Ef = 0.7 and Q c = 0.1, at T/ui ~ 0.2 the suscep- 
tibility is enhanced by a 8-9 % of its value at T = 0. 
Although this increase is rather small, it is nevertheless 
an interesting feature since it is not related to any in- 
crease of the lattice constant due to thermal expansion 
(the calculations reported here in fact are done for con- 
stant volumes). Experimentally in fact, x at T ~ 300 K 
is found to be larger than its value at T ~ 25 K of 
% and ~ 40 % for K 3 C6o and RbsCeo, respectively! 
Moreover recent data suggest a following decrease of x i n 
K3C60 by increasing temperature.ES The initial increase 
of x with temperature is currently explained by a temper- 
ature enhancement of the density of states at the Fermi 
level due to the thermal expansion of the unit cell. Note 
however that a power law dependence of the density of 
state on the lattice constant is not sufficient to reproduce 
the observed increase of x, while a stronger dependence 
like an exponential law needs a quite small value of the 
Stoner enhancement, p(J — Nol)^ 1 ~ 1.3, to fit the ex- 
perimental data of X-tZI Such small Stoner enhancement 
is in contrast to recent Maate Carlo calculations which 
estimate (1 — NqI)^ 1 — 30 Although a detailed study 
of the effect of the thermal expansion on x is beyond the 



scope of this paper, the results of Fig. [?] suggest that 
the electron-phonon contribution could be an additional 
source for the temperature dependence of x- Of course, 
a firmer evidence of the role of nonadiabaticity would 
be the measurement of the temperature dependence of x 
for a constant sample volume, in the spirit therefore of 
the constant volume resistivity experiment of some years 
ago.E] 

An alternative explanation of the non monotone tem- 
perature dependence of x has receatly been attributed to 
possible effects of mobile ions K + £3 This argument would 
predict however a finite isotope effect on x by alkali iso- 
topic substitution and no carbon isotope effect whereas 
in nonadiabatic theory an opposite trend is expected. 



V. CONCLUSIONS 

In this paper the nonadiabatic theory of the Pauli spin 
susceptibility in narrow-band systems has been formu- 
lated in order to identify possible signatures of nonadi- 
abatic electron-phonon coupling in fullcrides. We have 
identified peculiar features that can be considered as hall- 
marks of a relevant electron-phonon coupling in nonadi- 
abatic regime. In particular, an effective reduction of the 
spin susceptibility by the nonadiabatic electron-phonon 
coupling has been found in constrast to the conventional 
ME framework where no electron-phonon renormaliza- 
tion is expected. In addition, we predict a finite negative 
isotope effect on x which we suggest as possible experi- 
mental test. The Pauli spin susceptibility of x also ac- 
quires in nonadiabatic regime an anomalous temperature 
dependence in qualitative agreement with the experimen- 
tal data. 

The investigation of nonadiabatic effects in Cgo com- 
pounds attracts new and refreshed interest in the light of 
the recent experimental indications of an unconventional 
phonon-based origin of the superconducting pairing.Q 
The failure of the conventional theory of superconduc- 
tivity should not be surprising by considering that one 
of the basic assumptions of Migdal-Eliashberg theory, 
namely the adiabatic principle on which Migdal's the- 
orem relies, breaks down in fullerene compounds where 
phonon frequencies are of the same order of the Fermi 
energy. 

Aimed by these motivations, in the past years we have 
developed the nonadiabatic theory of superconductivity 
in narrow band systems that we propose as the theoret- 
ical framework to properly describe fullcrides and prob- 
ably cuprates. The inclusion of the nonadiabatic vertex 
corrections arising from the breakdown of Migdal's the- 
orem has permitted to explain in a natural way some of 
the experimental features of Cgo compounds which ap- 
peared anomalous within the ME context. For example, 
the low density of charge carriers can be regarded as a 
characteristic element since it is a direct by-product of 
the band narrowness and as a consequence of the nonadi- 
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abatic regime. In addition, the strong electron-electron 
interaction was shown to yield more attractive pairing 
channels than repulsive ones in nonadiabatic regime, and 
to play therefore a positive role with respect to the super- 
conducting onset. More technical anomalies, as the re- 
duction of T c in fullerides by induced disorder, have also 
received a natural explanation in the context of nonadi- 
abatic theory. 

The results of the present work shed new light also on 
the anomalous T c vs. x dependence in ammoniated alkali 
doped fullerides. In (NH^xA^Cqo indeed a sharp reduc- 
tion of T c , together with a weak expansion jaf the lattice 
constant a, is observed upon increasing of x.ErEl In addi- 
tion measurements of Pauli spin susceptibility show that 
such ra,|¥cduction of T c corresponds to an enhancement 
of This experimental situation looks quite puz- 

zling from the conventional point of view. In ME theory 
indeed an increase of the lattice constant a is reflected 
in an enhancement of the bare density of states at the 
Fermi level Nq. Both the Pauli susceptibility \ tx iVo 
and the critical temperature logT c oc —1/Nq are there- 
fore correspondingly expected to increase. Taking into 
account Stoner enhancement would make even stronger 
this trend. A theoretical explanation of this situation is 
still a open issue. 

The analysis of the Pauli susceptibility here presented 
suggests however a new possible interpretation. A cru- 
cial point is the observation that the Pauli susceptibility 
in the nonadiabatic theory, apart from the Stoner factor, 
does not give a direct probe of the electron density of 
states but it has to be considered as a quantity renor- 
malized by electron-phonon effects. Even more impor- 
tant is the fact that the renormalization itself depends 
on electron-phonon properties, so that an increase or 
decrease of x can be tune by electron-phonon proper- 
ties. Namely, in the model considered in this paper, the 
Pauli spin susceptibility x acquires a significant depen- 
dence on A and Q c . Figure |^ shows an enhancement of 
X by decreasing A or Q c (in vertex corrected theory). 
From an intuitive point of view such an enhancement 
of x should be accompanied by a depletion of T c . We 
quantify this concept by calculating the superconducting 
critical temperatures in non-crossing approximation and 
in vertex corrected theory. We consider here A = 0.7, 
Q c = 0.1, uj^jEp = 0.7 and N Q I = 0.4, and we study 
the T c vs. x dependence as function of A (left panel) and 
Qc (right panel). Numerical calculation are shown in in 
Fig. H for both the theories: solid lines represent the ver- 
tex corrected one and the dashed line the non-crossing 
approximation which does not depend on Q c . The re- 
sults shown in Fig. || suggest therefore that the driv- 
ing element of the anomalous T c vs. x behaviour could 
be the electron-phonon coupling in nonadiabatic regime 
and the weak lattice expansion of secondary importance. 
A microscopic relation between ammonia intercalation 
and possible modifications of the electron-phonon inter- 
action is still missing. The high critical temperature of 
(NH 3 )4Na2CsC6o however seems indicate that an impor- 
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tant role can be played by the displacement of the alkali 
cations from the center of the interstitial sites and by the 
splitting of the crystal field. 



APPENDIX A: 

Aim of this appendix is to provide the analyti- 
cal evaluation of the electron-phonon vertex correction 
P(k, k';n,m) and cross function C(k, k';n,m) and of 
their average over the momentum transfer. The follow- 
ing derivation is an improvement of the previous calcu- 
lations presented in the appendix of Ref. || by explicitely 
including the self-consistent self-energy renormalization. 
In order to do that a different procedure is employed. Let 
us consider first the vertex correction, Eq. (p6[), reported 
below for convenience: 

P(k, k'; n,m) = -T^2Yl - ^ D i n ~ l ) 

l q 

(Al) 



iWi - e(q) iWi- n+m - e(q - k + k') ' 

We consider a model for the electron-phonon matrix ele- 
ment which simulates the effects of the strong electronic 
correlation. Strong electron correlation has been shown 
in literature to favour the forward eLeetpan-phonon scat- 
tering at small transferred moment a. ll^Tta Following Ref. 
^, we therefore schematize the electron-phonon matrix 
| 5 (q)| 2 as: 



ls(q) 



MQc-Q), 



(A2) 



where Q = q/2kp, Q c = q c /2kp and kp is the Fermi 
wave vector. is the Heaviside step function and q c is a 
momentum cut-off which depends on the degree of cor- 
relation in the system: q c is smaller for more correlated 
materials. 

The explicit inclusion of the self-energy in W n does 
not allow for an analytical integration on the Matsubara 
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frequencies. In contrast with Ref. || we use thus the al- 
ternative procedure to evaluate first the integration over 
the momenta leaving the summation over the frequen- 
cies untouched. In evaluating the integration over the 
Brillouin zone, we employ some basilar approximations 
whic h are justified by the dimensionless cut-off Q c of Eq. 
(A2) supposed to be small (strong correlation case). In 
this context, we expand the electronic dispersion term 
e(p — k + k') as: 



e(p - k + k') ~ e(p) + 2E F Qa cos < 



(A3) 



where Ef = v F kF and a and <f> are the polar angles 
between q and k. In writing Eq. ( A3 ) , we have assumed 
a isotropic system and set |k| = jk 7 ] = |q| = kp, Q = 
|k — \s!\/2Kp. The smallness of the second term of Eq. 
( |A3| ) is enforced by the 0-functions of g(q) 2 present in 
Eq. ( |Al| ) and in front of the vertex function [see Eq. 

©]■ . . 

In similar wa y, th e leading order of the O-function ap- 
pearing in Eq. (Al) reads: 



0(<fc-|k-q|) ~6(2Q c -a). 



(A4) 



Finally, the sum over the Brillouin zone can be expressed 
in terms of the same angular coordinates a, (f>: 



E 

p 



2tt 



* da sin a 



N{e)de. 



(A5) 



-Eb 



By making use of Eqs. flAJ), <\A§, (JA5|) in Eq. ([Al]) 
we obtain the expression for the electron-phonon vertex 
function: 



P(k,k';n,m) = ^TV£>( 

Ef , l l 

de 

_ Ew iWi - e iW x 

—n+m 



2Q C 



— e — 2EpQa cos < 



daa 



(A6) 



where, coherently with the previous model, we have as- 
sume a constant DOS N(e) — No. 

It is convenient rewrite Eq. ( |A6| ) in the form: 

P(k,k';n,m) = A T ^D(n-Z) fj^ 

I 



2Qc daa 2E F Qa cos - i [Wi - Wi_ n+m ] 



2 [2P F Qacos</> 
e + iWi e • 



' + [Wi - W,_„ +m ] 2 
2E F Qa cos + iWi- n+m 



e 2 + Wf {e + 2E F Qacosi 



regimes Q/ [W; — Wi ] < I or Ql [Wi - Wi 
1 make arise the difference between the static and dy- 
namic limits and characterizes the main feature of the 
momentum-frequency structure of P. In contrast to that, 
the integrand in the third line of Eq. ( A7) is basically reg- 
ular in the whole momentum-frequency space. Hence, in 
the small Q expansion justified by the correlation cut-off 
parameter Q c , this latter term can be safely expanded at 
the second order in Q without losing important features 
of the vertex function. 

By integrating on e and expanding in Q we are then 
left with: 



P(k,k';n,m) 



2Qa daa A(n, m, I) - B(n, m, I) (2E F Qa cos < 



2 [2E F Qa cos ( 



[Wi - Wi- 



(A8) 



n+mj 



where 



A(n, m, I) 



{Wi - Wi^ n+m ) 
E F 



arctan 



arctan 



■n+m 



B(n,m,l) = (Wi-Wi- n+m ) 



E F Wi- 



Ef 
Wi 



■n+m 



(A9) 



[El 



W? I 2 

"i-n+mJ 



E 2 

F 



w? 



(A10) 



The last two integrals on a and </> can now be performed 
analytically giving: 

P(k, k'; n, m) = AP(Q, Q c ; n, m) 

= -XT D(n - I) {2B{n, m, I) 
i 



A(n, m, I) - B(n, m, I) [Wi - W X - 



n+mj 



(2E F QQ C 




2E F QQc 



Wi. 



n+m 



- 1 



(All) 



The average over the momentum Q is finally obtained 
from Eq. (p7|). Also in this case, the integration can 
be done analytically and the resulting averaged vertex 
function P(Q c ;n,m) becomes: 



P(Q C ; n, m) = Tj2D(n-l){B(n,m,l) 



(A7) 



A(n, m, I) - B(n, m, I) [Wi - Wi 



n+mj 



A very important role on the structure of the 
electron-phonon vertex function is played by the term 
in the second line of Eq. (A7). It contains a fac- 
tor 1/ [{2E F Qa cos <j)) 2 + (W t - W ; _„ +m ) 2 ] which is ev- 
idently nonanalitic for Q — > 0, Wi — ► Wi- n + m . The dif- 
ferent behaviour of this term with respect to the opposite 




(A12) 







The calculation of the cross function, Eq. (|T]), fol- 
lows basically the same lines reported above. The main 
difference with respect to the vertex function lies in the 
momentum dependence, which we handle in the same 
way as reported in Ref. ||[ Therefore, in Eq. ( |3l| ) we set 
g(k — q) 2 g(q — k') 2 ~ g(k — q) 2 g(k — k') 2 , which is an 
approximation valid for small values of Q c . Concerning 
the electronic dispersion, we expand e(k + k' — q) in the 
following way: 



e(k + k' - q) ~ e(q) + E F (1 - Q 2 )a 2 



-2E F Qy/l - Q 2 acos( 



(A13) 



where the angles a and 4> have the same meaning as be- 
fore. The integration over the energy and the angles can 
now be performed and, up to the leading order in the 
momentum transfer, the cross function becomes inde- 
pendent of Q = |k — k'l/2/ci?. Hence, for small values 
of Qc, the momentum average leaves the cross function 
unaltered and, from Eqs. (|32| , |33| ) , we obtain 

C(Q C ; n, m) = T ^ D(n - l)D(l - m) {2B(n, -m, I) 
i 

, A(n, -m, I) - B(n, -to, 1){Wi - W t -„- m ) 2 



x arctan 



2E F Ql\W l - Wi- n 
AE F Q 2 C 

\Wl - Wl-n-j, 



(A14) 



where th e functions A and B are reported in Eqs. 

(|a9|.[aTo|). 
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